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Abstract: We reconsider the strength of the electroweak phase transition in the sin- 
glet Majoron extension of the Standard Model, with a low (< TeV) scale of the singlet 
VEV. A strongly first order phase transition, of interest for electroweak baryogenesis, 
is found in sizeable regions of the parameter space, especially when the cross-coupling 
A/is|>S'p|-ffp between the singlet and the doublet Higgs is significant. Large Majorana 
Yukawa couplings of the singlet neutrinos, yiSu^Ui, are also important for strength- 
ening the transition. We incorporate the LEP and Tevatron constraints on the Higgs 
masses, and electroweak precision constraints, in our search for allowed parameters; 
successful examples include singlet masses ranging from 5 GeV to several TeV. Mod- 
els with a strong phase transition typically predict a nonstandard Higgs with mass 
in the range 113 GeV < niH ^ 200 GeV and production cross sections reduced by 
mixing with the singlet, with cos^ 6 significantly less than 1. We also find examples 
where the singlet is light and the decay H SS can modify the Higgs branching 
ratios relative to Standard Model expectations. 



*current address: Barclays Capital Japan Limited 



Contents 



|T|. Introduction 1 
g. The potential 3 
1^. Search of parameter space 6 



0| Algorithm 8 



Criteria for accepted points 9 



3.3 Distributions of parameters 12 



3^1 Regime of smaller couplings 15 



375| Constraints from electroweak precision observables 18 

1^. Implications for LHC 20 

m. Conclusions 25 

1^ . Analytic approximation for thermal potential 27 

|B|. Field dependent masses 28 

|C|. Renormalization constants 30 

|D]. Beta functions 30 

1^. Electroweak precision observables 31 



1. Introduction 

Nonstandard Higgs sectors are interesting from the perspective of LHC physics and 
cosmology. While the Standard Model predicts a smooth cross-over for the elec- 
troweak phase transition (EWPT) [1], extensions can give a strongly first-order phase 
transition, which is a necessary ingredient for electroweak baryogenesis, and could 
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also possibly generate observably large gravitational waves [2] or primordial mag- 
netic fields [3]. Supersymmctric extensions of the Standard Model have been studied 
the most intensively in this respect [4], but it is also possible to get a strong transi- 
tion from more generic two-Higgs doublet models [5,6], from technicolor theories [7], 
higher-dimension operators involving the Standard Model Higgs [8] , or from singlets 
which mix with the Standard Model Higgs [9-17]. In the last category, the singlet 
Majoron model [18] is an interesting example since it was originally motivated by the 
spontaneously breaking of lepton symmetry and consequent generation of neutrino 
masses by the seesaw mechanism. It is the model which we consider in the present 
work. 

In the singlet Majoron model, right-handed neutrinos v^^i acquire Majorana 
masses Mj = yi{S) through their Yukawa couphngs to the complex singlet field 
S, when it gets a VEV. Denoting the Yukawa couphngs to the doublet Higgs by 
hiUr^^iHun^i, the seesaw masses of the light neutrinos are given by 

where v is the VEV of H. If the Yukawa couphngs are 0(1), then I^S) ~ jmjj^i^ a 
very high scale > 10^^ GeV. This is the usual assumption, which would render the 
singlet field irrelevant for physics at the electroweak scale. However, we know that 
small Yukawa couphngs exist even in the Standard Model: that of the electron is 
O(10~^). If the hi are also of this order (while |/i ~ 1), then {S) could be as small as 
10~^^v^/(0.1 eV) ~ 300 GeV. Prom this point of view, a low scale for the singlet is 
no less natural than the Standard Model itself, and merits consideration. 

The effect of the singlet Majoron on the EWPT has been considered previously 
in [17], but these papers were written before the final LEP/Tevatron bounds on the 
Higgs boson mass van oi' values of electroweak precision obscrvables (EWPO) were 
known, and thus they could not take these important constraints into account.^ As 
is well appreciated, the strength of the EWPT tends to be inversely related to van] 
moreover the EWPO constraints tend to exclude heavy singlet fields which have 
significant mixing with the Higgs doublet. There have also been studies of related 
models [10-16], where the singlet is a real field, or a complex one such that the global 
U(l) symmetry under which S might transform is explicitly broken by terms like . 

^It is possible to evade the EWPO constraints, so we will present results both with and without 
applying them. 
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These models are also very interesting, but sufficiently different from the Majoron 
model to justify a separate study of the latter. The models with explicitly broken 
symmetry are more generic, but not motivated by considerations of neutrino physics. 
Moreover, cubic terms in the scalar potential tend to make it easier to find a ffist- 
order phase transition, so we would expect the physics leading to a strong EWPT 
to be qualitatively different in the two classes of models. Indeed we will show that 
the coupling of the right-handed neutrinos to S plays an important role in getting a 
strong phase transition in the Majoron model. 

Because the global lepton symmetry is spontaneously broken, the imaginary part 
of S* is a Goldstone boson, the Majoron. Since we are assuming the doublet Yukawa 
couphngs to be quite small {hi < 0(10^^)), the massless Majoron couples very weakly 
to the light neutrinos, with strength hiv/{S). These couplings are diagonal in the 
mass basis at this order; off-diagonal couplings which could lead to neutrino decays 
are suppressed by {hiv/{S)Y (see ref. [19] and references therein). Such a weakly 
coupled Majoron goes out of equilibrium well before nucleosynthesis, and also has a 
negligible effect on energy loss from stars, and so it is experimentally unconstrained.^ 

In the remainder of the paper, we derive the finite-temperature effective potential 
(section and present our methods and results for the strength of the EWPT from a 
wide search of the model's parameter space (section H). We analyze the nonstandard 
decay modes and discovery potential of the singlet sector at the LHC in section 
^. Conclusions are given in section |[ We present formulas for field-dependent and 
thermal masses needed for the potential in the appendices, as well as those pertaining 
to our renormalization prescription, the running of the couplings, and formulas for 
the oblique parameters for electroweak precision observables. 



2. The potential 

Because the left-handed neutrino Yukawa couplings are assumed to be very small 
for our purposes, we neglect them in what follows. Similarly, only the top quark 
is retained amongst the other fermions of the Standard Model. At tree level the 

^In contrast to the scenarios discussed in ref. [15], we do not require the singlet to provide a 
dark matter candidate. 
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c. 



potential is then 

(2.1) 

in terms of the complex Higgs doublet H = [H^, H^), complex singlet S, top quark 
and right-handed neutrinos. For definiteness, we take three generations of right- 
handed neutrinos with equal Majorana Yukawa couplings yi. Due to the cross- 
coupling A/js, Vh and Vg are not generally the VEV's of the fields at the minimum of 
the potential. Rather, the relation is 

vl = 2{Hr + ^{Sf; (2.2) 
Aft 

We take (H) = 174 GeV and (S) (to be varied) as the physical input parameters. 
Because of the Z2 symmetries H —H and S —S in the scalar potential, there 
is no loss in generality in assuming that both VEV's are positive. (The signs of the 
fermion masses are not physically significant.) 

At finite temperature, the lowest order thermal correction to Vq is a function of 
the field- dependent particle masses, mi{H, S): 

( bosons 
\ fermions 



AVT = Ty± [ /4ln f 1 ^ e-Vp'+'"^(^'^) 
4- J (27r)3 V ^ 



where (3 = 1/T. These functions are often approximated by their high temperature 
expansions, but for numerical purposes it is preferable to use an approximation that 
works at all values of rrii/T. We use the approximation described in ref. [5], in which 
the high-T and low-T expansions are smoothly joined together at some large value 
of rrii/T. This is reviewed in appendix The expressions for the field-dependent 
thermal masses are given in appendix 0. 

Furthermore, it is important to improve the thermal contribution by resumming 
the ring diagrams, which amounts to replacing with the thermally corrected 
masses (of the form + qT^) in eqs. ( |2.3| ). Otherwise there is a danger of over- 
estimating the strength of a first order phase transition. Often, this substitution is 
made only in the cubic term of the high-T expansion, where it has the biggest effect 
on the barrier between the true and false vacua of the potential [20]. However, when 
one uses an expression that correctly captures both the large and the small rrii/T 
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behavior of (|2.3|) , there is no way to consistently include thermal mass effects only 
in the cubic term, since it appears explicitly only in the high-T expansion. Trying 
to do so creates a kink in the potential when the high-T and low-T expansions are 
joined onto each other. To avoid such complications, we simply replace mf by its 
thermally corrected value everywhere in ( p.3|) [21]. The discrepancy between the two 
approaches is formally significant only when computing the potential to two loops. 

The thermal potential ( |2.3| ) can be regarded as a one- loop effect, so for con- 
sistency one must also include the one-loop, zero-temperature (Coleman- Weinberg) 
correction to the potential, 

ATA l^lrrl2 1 4/rr r.N A ^liH ■, S) 3\ f+1, boSOUS 

AVc„. = -A\H? + — Y^ :,4(E, S) \X. - 2 j M -1, tennio„= '^■^> 

where Ajifp is a counterterm and /i is the renormalization scale. We do not include a 
counterterm of the form iJiS'p because this can be absorbed into a redefinition of v^- 
However for the doublet Higgs it is convenient to introduce the Ajifp counterterm, 
because then one can maintain the tree-level relation between the Higgs VEV 
and Vh- In fact, it is convenient to maintain both relations (|2.2| ), so that the position 
of the zero-temperature minimum of the potential H = (H), S = (S), is known 
analytically. We thus adopt as our renormalization prescription 

dV dV 

— = — = at H={H),S={S) (2.5) 

where V = Vq + AVcw- These two equations can be analytically solved to find /i 
and A (see appendix 0). 

A notable feature of AVcw is that the terms \n{mf{H, S)) appear in such as way 
as to exactly cancel corresponding terms in the high-T expansion of AVr] then the 
cubic term (m^)^/^ of AVr is the only source of nonanalytic dependence on the fields. 
To preserve this property, we also replace mf by the thermally corrected expression 
in AVcw when we do the ring improvement of the potential [21]. 

A complication which arises in the effective potential is the appearance of neg- 
ative values of mf{H, S) for the Goldstone boson degrees of freedom at small H or 
5"; this can happen even when the thermal correction to is included. Such values 
create a problem with the cubic term (m^)^/^ in the high-T expansion of the ther- 
mal potential. Even if one takes only the real part, which vanishes for negative m^, 
derivatives of this with respect to the fields are discontinuous at the point where 
changes sign, leading to serious difficulties for algorithms which attempt to minimize 
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the potential. There are various prescriptions in the hterature for deahng with the 
Goldstone bosons. We take the simplest approach, which is to simply omit their 
contributions from ISVcw and ISVt- Experience with other models indicates that 
the Goldstone bosons never have a strong effect on the phase transition in any case. 
Thus we omit the contributions from if^, Im(if°), and the Majoron j in the sums. 

In the limit of very heavy singlet and neutrinos, one expects the effects of the new 
physics to decouple. This is evident in the low-T expansion of the thermal correc- 
tion to the potential, since the effects of heavy particles are Boltzmann suppressed. 
However, this decoupling is missing from the naive thermal corrections to the masses, 
( P.4| - |B.4D , which were derived from the high-T expansion of the potential. To correct 
for this, we insert Boltzmann factors involving the heavy particle masses, 

bm\ = Qa. + ^A..e-«/^ + 1(3/ + g'^) + \y^^ (2.6) 
Sml = (^IXs + lx,s + 1 E y^-™^/^ j (2.7) 

where rris and are evaluated at zero temperature and at the minimum of the 
potential. This procedure is somewhat rough, but better than ignoring the issue 
altogether [22]. 



3. Search of parameter space 

Our goal was to make a broad scan of the parameter space, in search of models giving 
a sufficiently strong EWPT for electroweak baryogenesis. We used a 30 x 15^ x 20 x 30 
grid (~ 4 X 10^ points) on the five parameters Xhs, Xh, Xg, Ui (taking equal Majorana 
Yukawa couplings for 3 generations of right-handed neutrinos) and (s), the VEV of 
the real component of S*, (s) = ^/2{S). After a preliminary scan of the parameter 
space to determine the values of interest, these were taken to be in the ranges 

< A„ Xh < 3, -3 < Xhs < 3, 0<yi < 8, (s) < 1800 GeV (3.1) 

subject to the constraint Xhs > — V A^A^ which is needed for V to be bounded from 
below for large field values (at tree level). The chosen ranges include large values of 
the coupling constants, but since Xf / Att , yf / An < 1, they are not unreasonably large. 

We find that the typical pattern of symmetry breaking is to first develop a finite- 
temperature VEV at 5 = Sci for S alone as T is lowered from very high values, i.e., 
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Figure 1: Trajectories of typical EWPT in H-S field space. The significance of the mixing 
angle 6 is discussed near eq. (3.9). 

the phase transition for the singlet to condense usually occurs above the EWPT. 
The electroweak transition is a jump from this false vacuum along the S axis to a 
true vacuum in which H ^ 0; furthermore (S) = Sc2 increases relative to its value 
in the (H) = minimum if X^g < 0, and decreases if Xhs > 0. This pattern is shown 
in figure 0. Although Sd is typically large, we will see that it can sometimes (when 
Xhs < 0) be zero. 

When \Xhs\ < 1, the mixing between the H and S fields is small, and the transi- 
tion takes place mostly in the H direction. For larger values of \Xhs\, which are more 
typical of cases with a strong EWPT, the induced H-S mixing causes the tunneling 
path to be along a linear combination of the fields. This is illustrated in figure ^ 
It is easy to understand the shapes of the valleys in the potential indicated in figs. 

First, the discrete symmetry H —>■ —H implies that the light direction will be 
purely parallel to the H axis at if = 0. When H gets a VEV, the term A/js|-f^P|5'p 
makes the squared mass of S more negative if Xhs < 0, increasing the VEV of S in 
the electroweak symmetry breaking vacuum. Conversely the VEV of S is decreased 
for Xhs > 0. 

The above discussion makes it possible to understand why the mixing term 
A/is|iJp|S'p can generally strengthen the EWPT. Roughly, we expect the VEV of H 
to scale like the square root of the negative mass squared term for H. When S 
has a VEV, — /i^ gets a contribution XhsiS)"^. Now suppose that (5*)^ changes by the 
amount 6{S)^ when H makes the transition between the symmetric and electroweak 
symmetry breaking EWSB vacua. We expect that the critical value increases 
with Xhs S{S)'^. The preceding discussion shows that this quantity is always positive. 
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Figure 2: Shape of potential V{H, S) at the critical temperature for large value of Xhs = 
—2.4 (left) and small value X^s = —0.2 (right), illustrating the effect of Xhs on the path in 
field space between the degenerate minima. Contours of the potential are projected onto 
the lower plane. H, S are in units of 100 GeV. 

regardless of the sign of Xhs, so the mixing term should tend to strengthen the EWPT 
whenever it has a significant size. 

3.1 Algorithm 

Because the phase transition typically proceeds in two steps, due to the two con- 
densing fields, automating the search for a strong first order transition proved to be 
somewhat more difficult in this model than for effectively single-field models. The 
key steps are to identify whether there is a barrier between the trivial (H) = min- 
imum and the EWSB (H) ^ minimum, and to bracket the critical temperature if 
there is one. Figure ^ outlines our algorithm. It outputs a logical variable success 
to indicate whether a first order transition with Vc/Tc > 1 was found, and bracketing 
temperatures Tmin and Tmax if so. 

The algorithm shown in figure ^ is simple to implement in a single-field model, 
but when there are two fields it can be difficult to properly identify the relevant 
direction along which to check the curvature of the potential. In the present case, we 
start by finding the global electroweak symmetry breaking minimum near T = 0, and 
a local S-breaking minimum or saddle point on the 5* axis, and raise the temperature 
until these two critical points become degenerate. Naively, the tunneling path would 
be along a line connecting these two points, but because of the "banana" effect 
shown in fig. ^ the relevant flat direction might really be curved. One must therefore 
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Algorithm to bracket critical temperature of first order phase 
transition between How and Thigh, success = .false, if there is 
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Figure 3: Flowchart for algorithm to bracket the critical temperature for a first order 
phase transition. 

determine the curvature locally around the two putative minima to check that they 
really are minima. If the transition is actually second order but the transition path 
is curved, one could mistakenly conclude that it is first order by measuring the 
curvature of the potential along the straight line rather than the curved path. Of 
course a visual inspection of the shape of the potential would eliminate such cases, 
but we needed to automate this. To do so, we minimize the potential on small circles 
surrounding the putative local minima, to verify that they indeed are minima, and 
to find the directions of shallowest ascent. 

3.2 Criteria for accepted points 

The basic requirement for a strong enough phase transition is 

^ ^ 1 (3.2) 

c 
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Figure 4: 95% c.l. LEP bound on mixing angle (cos^ 6) of the mostly-doublet Higgs state, 
from table 14 of ref. [25] . The same bound applies to sin^ 9 for the mostly-singlet state. 

where Vc is the VEV of the real Higgs field h {^/2 x 174 = 246 GeV at zero temper- 
ature) at the critical temperature Tc [23] (see [24] for a pedagogical review). This 
avoids the washout of baryons produced during the EWPT by sphalerons. A more 
careful treatment would be to calculate the sphaleron energy in the model at hand, 
since this can in principle be different from the Standard Model value and change 
the bound. The change is typically small however, and so we do not consider this 
effect. 

In addition, we demand that the LEP limit on the Higgs mass [25] be satisfied. 
In this regard, another important feature of the A/is|i/^| 15*^1 interaction is that it 
can cause large mixing between the singlet and doublet Higgs bosons, leading to a 
reduction in the production cross section for the hghtest mass eigenstate [26]. We 
take the fluctuations of the flavor and mass eigenstates to be related via 

ss J \-sm9cos9 J \ SS' J 

We restrict cos 9 > 1/ \/2, so that H' is the "Higgs-like" state and S' is the "singlet- 
like" state, regardless of which one is heavier. The mixing angle suppresses the 
couphngs of either state relative to the couphngs of a SM Higgs boson. The pro- 
duction cross section of the Higgs-hke state is reduced by cos^ 9, while that of the 
singlet-like state scales hke sin^ 9. We demand that both of these are less than the 
LEP limit; i.e., both cos^^ and sin^^, evaluated at the appropriate mass, must be 
less than the value in column (a) of table 14 of ref. [25]. The bound is shown in 
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figure |. The procedure of applying tlie same bound independently to botli states 
must be modified if tlie two are close to each other in mass. However we will find 
that for cases that give rise to a strong EWPT, there is always a large separation 
between the masses, justifying this simpler approach.'^ 

Recently, the DO and CDF experiments have disclosed new limits excluding the 
SM Higgs boson in the region 160 — 170 GeV [27]. We fit the limit on cos^(^) 
from Table XIX in this region with the quadratic function cos^ ^max = 142.43 — 
1.716 m//' + 0.0052 m|^, where tuh' is in GeV. This has values 0.99, 0.86, 0.99 at 
rriH' = 160, 165, 170. Because of this limited range, the effect is small in our first 
search of parameter space, which involves a large range of couplings. However in the 
search which targets smaller couplings, the CDF/DO limit covers a larger fraction 
of the range of allowed masses, and the mixing angles tend to be smaller, and so 
the constraint has a more pronounced effect: 450 out of 5300 parameter sets are 
removed. 

Some cases of interest have very light singlets. Ref. [28] noted that values of 
rris' < 5 GeV (the B meson mass) are strongly constrained {6 < 10^^) by the decays 
B S'X followed by S' /^^A*^- We thus exclude nis' < 5 GeV if the mixing 
angle is greater than 0.01. This has a negligible effect on our broader search of the 
parameter space since a very small fraction of this sample has light singlets, but in 
the search which is limited to smaller values of the coupling constants, this constraint 
is more significant. 

Another important criterion is that the EWSB vacuum at T = must exist. 
Although this seems obvious, our broad scan of parameter space includes cases where, 
due to radiative corrections to the tree-level Higgs potential, the curvature of the 
potential is positive when H = 0, leading to no EWSB. In fact, we will see that some 
concentrations of the preferred parameter space tend to be close to this perilous edge, 
especially when X^s < 0. 

Since we consider models with large couplings and large masses, a consistency 
requirement for perturbation theory to be under control is that none of the running 
couplings diverge (reaching a Landau pole) at renormalization scales smaller than 
the heaviest particle masses. The beta functions for the largest couplings are given 
in appendix 0. For each otherwise accepted parameter set, we integrate these to find 

^In particular, we find no cases where both niH' and ms' < 114.4 GeV; in this situation the 
decays H' S'S' or S" H'H' (if kincmatically allowed) could modify the branching ratios with 
respect to the Standard Model prediction. 
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the first Landau pole and discard parameters which fail this test. This eliminates 
approximately 15% of otherwise accepted parameter sets from the range ( |3.1| ). 

Finally, we include constraints on the oblique parameters S, T, U from elec- 
troweak precision observables (EWPO). In order to not to mask the intrinsic dy- 
namics of the phase transition too much, we chose to first present results without 



inclusion of the EWPO constraint. A separate section p.5| is devoted to showing how 
the results are affected by its inclusion. We give details about its implementation 
there and in appendix |^. 



3.3 Distributions of parameters 

Out of the 4 X 10^ points tested on our uniform grid in the parameter space, approx- 
imately 0.07% generate a strong enough phase transition and fulfill the other criteria 
mentioned above. We display the distributions of accepted parameters in figure ^ 




0.5 0.6 0.7 0.8 0.9 i 5 10 15 1 1-5 2 2.5 



Figure 5: Distributions of parameters which generate a strong enough first order phase 
transition. Lighter (red) bars correspond to X^g < 0, darker (black) to X^g > 0, and \Xhs\ 
is shown in the distribution for X^s. it sign indicates regions associated with X^g ^0. {s) 
is the real (not complex) singlet VEV. 
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The samples are divided into two groups, according to whether Xhs < or Xhs > 0, 
due to the expected quahtative differences between the two cases. These differences 
are highhghted by the separate distributions shown for Xhs < and Xhs > 0. One 
feature which they have in common however is the need for generally large values of 
\Xhs\, in agreement with our argument that H-S mixing is important for boosting 
the strength of the phase transition. 

Another striking feature of the distributions is the preference for large values 
of the Majorana neutrino Yukawa coupling, yf, and the largest values being corre- 
lated with Xhs > 0. It was pointed out in ref. [29] that new heavy fermions with a 
large Yukawa coupling to the Higgs could strengthen the EWPT. One might wonder 
whether this could be the origin of the need for large yf in our model, since both 
VEV's S and H are changing during the phase transition. However, the sign goes 
the wrong way. The heavy fermion effect requires the fermion to be heavier in the 
EWSB phase than in the symmetric phase. Fig. ^ shows that this is the case when 
Xhs < 0, but not for Xhs > 0. However the preference for large yf is greater for 
Xhs > in fig. |. 

To understand this behavior, we have varied yf away from the accepted value 
for some sample points. We typically find that Vc/T^ is an increasing function of yf, 
while the light Higgs mass is decreasing. There is thus a tension between the demand 
for a strong phase transition and the LEP bound, which results in a narrow window 
of allowed yf values, while keeping the other couplings fixed. These dependencies on 
yf are illustrated around a sample accepted point in fig. §(a). 

By testing many hypotheses, we eventually discovered an analytic explanation for 
the trends visible in the distributions of parameters in fig. |^. It depends crucially on 
the one-loop zero-temperature correction to the effective potential. Let us try to make 
a rough analytic estimate of the strength of the phase transition, which is charac- 
terized by Vc/Tc. The critical temperature is approximately where the temperature- 
dependent mass squared of H, evaluated at {H,S) = (0,S'ci), goes through zero. 
Using the field-dependent mass m^^ of ( p.l|) , the one- loop counterterm A of ( |2.4| ), 
and the temperature correction ( [B.4| ), we get Tf ~ (m^^(0, Sd) — A) /{6mff^/T^). 
Putting these results together gives the estimate 



The critical temperature can become small relative to Vc if the Higgs mass renormal- 



VXh{H)^-Xhs{Sf,-{Sy)-A 



(3.4) 
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Figure 6: (a) Left: Tc, lightest Higgs mass rwH" and v^/Tc as a function of y^, for 
representative point A^^ = 1.90, = 1.70, A., = 0.90, (s) = 570 GeV. (b) Right: the r.h.s. 
of eq. (3.7) versus yf for parameter sets leading to a strong EWPT. Solid line is the fit to 
the cluster of points which nearly satisfy this linear relation. 

ization constant A, eq. ( p.2| ), becomes large. This can happen when the renormal- 
ization scale In/i^, eq. ( P.l|) , becomes large. It is straightforward to show that 

^ ^ (24A2 - 6yf + XI) {Sy + (AL + A..(4A. + 6A.)) (H)^ ^ ' > 

where the numerator is correct in order of magnitude, while the denominator is exact. 
Since In/i^ appears in A, eq. ( |C.2| ), large values of In/x^ can cause A to be large, 

"^-TK^ (^(^^' +^L)(^)' + A..O(A„ A., X,s){S)') (3.6) 
The denominator of eq. (|3.5|) vanishes when the relation 



yf = I - gK, + + 4AftA,..) I (3.7) 



is satisfied. The correlation between yf and the r.h.s. of eq. (|3.71) is shown in fig. 
P(b). This relation is just what one would expect from trying to minimize Tc by 
making In/x^, hence A, large. For smaller yf, the strength of the transition rapidly 
diminishes. For larger values, we lose the EWSB vacuum because the curvature mlf^ 
has the wrong sign at H = 0. 

This effect also allows us to understand other features of the distributions shown 
in fig. ^. Eq. (|3l7| ) shows that when Xhs > 0, larger values of yf result due to the 
term AXhXhs] this trend is seen in the histogram for yf. Similarly, this term puts a 
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limit on the magnitude of Xh when Xhs < but not when Xhs > 0, in agreement with 
the histogram for Xh- Moreover, inverting the relation (|3.7| ) to express in terms of 
the other variables readily explains why larger values of As are favored for A^s < 0. 
Solving for {H)"^ j {S)"^ similarly shows why larger values of {S) occur for Xhs > 0. 
In short, the relation ( p.7|) allows us to qualitatively understand most of the trends 
exhibited in fig. ^. The Coleman- Weinberg potential thus plays an important role 
in strengthening the EWPT. We have further tested this conclusion by running our 
program with the one-loop zero-temperature correction turned off, finding that the 
number of accepted points is drastically reduced in this case. 

- — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 1 We mentioned earlier that the 
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Figure 7: Distribution of Sd, the scalar VEV in 
the electroweak symmetric vacuum at the critical 
temperature. 



VEV of S in the electroweak sym- 
metric vacuum, ^d, can also be zero, 
but only when Xhs < 0. The dis- 
tribution of Sci in figure |^ indeed 
shows a spike at Sd = for Xhs < 0. 
Fig. 1^ makes it clear why there such 
is a correlation between Sd and Xhs'- 
for Xhs > 0, S must decrease during 
the transition. It cannot do so if it is 
already zero. Eq. ( |3.4| ) also gives in- 
sight into this correlation: large val- 
ues of Sd are disfavored when Xhs < 



since this tends to increase and decrease Vc/Tc. 



3.4 Regime of smaller couplings 

It is interesting to know how important it is to have large coupling constants to 
get a significant effect on the EWPT. With such large couplings as the maximum 
values used in the scan described above, one expects to reach a Landau pole in one 
of the couplings at a relatively low renormalization scale /Umax- We have investigated 
this by integrating the renormalization group equations for Xh, Xs, Xhs, Hi and yt 
(appendix |D|). The distributions of /imax plotted in figure || show that indeed new 
physics beyond the singlet and the right-handed neutrinos must typically come in 
at the scale of several TeV. (This would also be true in any model which uses the 
mechanism of ref. [29], since strong couplings are always needed to get that effect.) 
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This motivated us to further explore the model at somewhat weaker interaction 
strengths. We thus made an additional scan of the same point density as described 
previously, but in the more limited range 

< A„ A,, < 1, -1 < \hs < 1, 0<y^ < 2, (s) < 1800 GeV (3.8) 

Encouragingly, none of the accepted parameters from this range fail the Landau pole 
test, showing that indeed perturbation theory is more reliable in this case. (In fact 
/imax is always greater than 50 TeV for this sample.) We find that almost no points 
are accepted for X^s < 0, nor for (s) > 100 GeV, but there is a sizeable number 
of accepted parameters with Xhs > and (s) < 50 GeV. These points have small 
mixing angles, sin^^^ < 0.02 (due to the LEP constraint), and small singlet masses, 
ms ^ 20 GeV. Their phase transitions tend to follow a circular arc in the H-S plane, 
between the H and S shown in fig. |^(a). A curious feature is that Sd 

and Sc2 are typically an order of magnitude larger than the small zero-temperature 
VEV (S), due to (H) and hence the mixing being different at high temperature 
relative to T = 0. The parameter distributions are shown in figure 0. Unlike in 
the broader region of parameter space described in the previous section, here we do 
not find any correlation like that in eq. ( |3.7| ); thus these points give rise to a strong 
phase transition for different reasons than the majority of those in the large coupling 
regime. 

The dependence of the EWPT on the Yukawa coupling yf is due to the strong 
influence of yf on the 5* dynamics, which subsequently affects the dynamics of H 
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Figure 8: Distributions of maximum renormalization scale, where Landau pole develops, 
for the large-coupling parameter sets (left panel) and the smaller-coupling ones (right 
panel). 
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Figure 9: (a) Left: potential at critical temperature for typical point having smaller 
couplings: Xhs = 0.633, Xh = 0.633, = 0.167, yf = 1.25, {s) = 10 GeV. (b) Right: 
Dependences of Tc, Vc, Vc/Tc, Sd, and Sc2 on around the point specified in (a). 

through mixing. To understand this, we examined the dependences of various quan- 
tities, Tc, Vc, Vc/Tc, Sci, and Sc2 upon the Majorana Yukawa coupling yf. An example 
is shown in figure |^(b). There is a notable rise in Vc associated with the decrease 
in Sc2 — recall that this is the value of (S) in the EWSB vacuum at T = T^, see fig. 
|l|. We can give an analytic explanation for the relation between Vc and Sc2- To this 
end, it is useful to think in terms of an effective potential along the light direction 
H', which we approximate by the straight line paths connecting the symmetric and 
EWSB vacua shown in figure 0. At the critical temperature, we can write 

(3.9) 





where the mixing angle 6 is generally different from the zero-temperature mixing 
angle 6. At T = Tc, the shape of the potential is roughly of the form X'H''^{H' —v'^y = 
XH'^ — 2g'H'^ + fi^H'"^. From this form, we see that v'^ = g'/X'. The cubic term can 
be estimated from the tree-level potential as 

V,^uc = -2sg {Xhs4 + 2Xssl) SciH'^ ^ g' = ss{Xhs4 + 2Xssl)Sci (3.10) 

Similarly, the effective quartic coupling is A' = c|A/i -|- 2c|s|A/js + s'^Xg. This gives the 
estimate 

(4Ah + 2cpiXhs + s^Xs) 
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Figure 10: Distributions of accepted parameters with smaller values of the coupling 
constants, and Xhs > 0. 



From eq. and fig. |l] it is clear that for X^s > 0, a decrease in Sc2 leads to 

an increase the mixing angle 6 and consequently an increase in Vc- This accounts for 
the initial growth in Vc for yf < 0.6 in fig. ||. Beyond this point, Sc2 remains constant, 
but Sci decreases (hence 9 decreases), leading to a decrease in Vc- At the same time, 
eq. ( p.4|) shows that, for Xhs > 0, decreasing Sd leads to an increase in Tc. Both of 
the these effects cause Vc/Tc to go down with yf as observed in fig. ^(b). 

3.5 Constraints from electroweak precision observables 

In refs. [14, 30] it was noted that electroweak precision observables provide a strong 
constraint on the related model containing a real singlet field. It is known that the 
oblique parameters S, T, U are best fit by a light Higgs boson, and this preference 
thus extends to singlets that mix with the doublet Higgs. The same constraints as for 
the real singlet apply to the Majoron model, since the extra Goldstone boson does 
not mix and therefore plays no role. We have thus carried out the same analysis 
as in [14, 30] to further constrain the accepted parameter sets described above. For 
completeness, the relevant formulas are given in appendix |^. As a check on our 
implementation, we reproduced the results shown in figures 9 and 10 of [14]. 
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The EWPO constraint indeed has a strong impact on the accepted parameter 
distributions. In our larger couphng sample, 730 out of 1000 points are removed for 
Xhs < 0, and 1650 out of 1710 are excluded for Xhs > 0; overall 88% of otherwise 
accepted points are thus ruled out, at 95% c.l. The resulting distributions are shown 
in fig. [Tl|. The most striking difference relative to the corresponding results without 
EWPO, fig. is the elimination of large values of the doublet-like mass ttlh', and the 
restriction to smaller mixing angles, cos^^^ > 0.8. There is also a stronger exclusion 
of small values of the Majorana Yukawa coupling. 

In the small-coupling sample, the EWPO constraint is even more powerful, elimi- 
nating 4370 out of 4860 parameter sets. Again, the effect is to eliminate higher values 



of the Higgs boson mass. The distributions are plotted in fig. |I2|. The lower niH' 
values lead to correspondingly smaller values of the coupling Xh. 

Two cautionary remarks are, however, in order here. First, since the additional 
singlet can be light, a precision analysis of EWPO would require the inclusion of 
further parameters U, V, X [31,32]. Here, we do not go beyond S,T,U but leave 
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Figure 11: Distributions of accepted parameters with larger values of the coupling con- 
stants, after EWPO cut. 
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Figure 12: Distributions of accepted parameters with smaller values of the coupling 
constants, after EWPO cut. 

an extensive EWPO analysis for future work. Second, the presence of additional 
new physics beyond the singlet-extension of the SM could considerably weaken the 
EWPO constraints. Therefore we refrain from imposing EWPO as a strict constraint 
in the following. Rather, we present results both with and without this constraint. 



4. Implications for LHC 

From the perspective of collider phenomenology, it is quite intriguing that all our 
accepted points feature a relatively light scalar, either the singlet- or the doublet-like 
state, with mass less than about 200 GeV. The other state is typically considerably 
heavier. Concerning detectability at the LHC, it is important to to know how the 
mass of the lighter state correlates with the mixing angle, i.e., how its couplings 



compare to those of a SM Higgs boson. To this purpose, figure |I3| shows the scale 
factor C, of the squared couplings versus mass of the lighter scalar. Specifically, 
^ = cos^ 9 for the Higgs-like boson, while ^ = sin^ 9 for the singlet-like state. This 
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Figure 13: Distribution of scale factor of the squared couplings (relative to a SM Higgs) 
versus mass of the lighter scalar. For a doublet-like state ^ = cos^ 9 > 0.5, while for a 
singlet-like state ^ = sin^ 9 < 0.5. The vertical red line indicates the limit on a SM Higgs 
boson. The left (right) panel is without (with) the EWPO constraint. 



can be compared with the LEP bound (fig. ^) to get a feeling for how easily detectable 
the light boson may be at the LHC 

Let us first discuss the situation without the EWPO constraint, shown in the left 
plot of figure The density of points in this scatter plot indicates that there would 
have been many examples providing a strong EWPT in the LEP-excluded region, 
but there is not a strong bias toward being close to the limit. Nevertheless, there 
is an upper limit of mm{mH', rris') ^ 200 GeV from the requirement of a strong 
EWPT (note also the lower limit on the mostly-doublet state of niH' ^ 113 GeV 
from LEP data). Moreover, a large fraction of the accepted points features a sizable 
doublet-singlet mixing. 

In the sample with small values of the couplings, the situation is different, because 
only small values of < 20 GeV, and the mixing angle sin^ 9 < 0.02, are present. 
This situation was considered for the similar model of a real singlet field in ref. [28]. 
There it was noted that values of < 5 GeV (the B meson mass) are strongly 
constrained {6 < 10^^) by the decays B S'X followed by S' fi'^^~ . For ms' in 
the range 5 GeV < rris' < 50 GeV, the Higgs can decay into singlets, H' — > S'S' at a 
level which can compete with the two photon final state, H' 77. For the accepted 
parameters in the small coupling regime, we find that the right-handed neutrinos 
are always heavier than ■ms>/2. Therefore there is never an invisible decay channel 
S' vrVr in this case. Instead, S' decays predominantly into hb quark pairs due to 
the small doublet-singlet mixing. 
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(left) and with (right) EWPO constraint. The full red lines (horizontal and vertical) 
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The situation changes quite drastically when applying the EWPO constraint, 
as shown in the right plot of figure In this case, the allowed range shrinks to 
mm{mH', rns') < 156 GeV, and the regions with large mixing and/or a light singlet 
are almost completely cut away. 

It is also useful to consider the correlation between the mostly-singlet and mostly- 
doublet Higgs masses, shown in figure |T^, which reveals that typically one state is 
significantly lighter than the other. For both, \hs < and \hs > 0, we see distinct 
islands with either (i) rriH' ^ 200 GeV and m^/ ranging from few hundred GeV up 
to order TeV, or (ii) rriH' > 200 GeV and a (much) lighter singlet. The latter region 
is, however, completely removed by the EWPO constraint. For Xhs < 0, there also 
exists a small region with rriH' ~ 115 — 200 GeV and ms' ~ 50 — 130 GeV. The triple 
correlation between the two masses and the mixing angle is shown in figs. |T^ and |TB[ 



Regarding the decay modes, we see from figure |1J that for most points Higgs-to- 
Higgs decays, either 5" H'H' or H' S'S', are allowed, which could modify the 
branching ratios relative to those of a SM Higgs boson. (If Higgs-to-Higgs decays are 
absent, the branching ratios of both the H' and the S' are just the same as those of 
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Figure 15: Triple correlation of mass eigenvalues and mixing angle for Xhs < (left) and 
Xhs > (right); without EWPO constraint. 
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Figure 16: Same as figure 15 but with EWPO constraint applied. 



a SM Higgs.) Denoting the lighter of the two states as hi and the heavier one as h2, 
the generic expression for the decay width is 



T{h2 hihi) -- 
where cubic coupling g2ii is given by 



9211 
8711712 



1 — Amf/ml , 



(4.1) 



gn'S'S' = 24:Xh{h)cgsl-24:X,{s)sQcl+4:\hs {{h)ce{cl - 2sl) - {s)se{sl - 2c])) (4.2) 

for h2 = H', hi = S' and 

gs'H'H' = 2AXh{h)sQcl + 2AXs{s)cesl + AXhs {{h)ss{sl - 2cl) + (s)ce(c^ - 2s^)) (4.3) 

for h2 = S', hi = H' . Here, cq = cos 6 and sq = sinO. The branching ratio is then 
given by 



BR{h2 hihi] 



r(/i2 — > hihi 



(4.4) 



^T{HsM) + Tih2^hihi)' 
where ^ = cos^ 9 (sin^ 9) for h2 = H' (S") and r(ifsM) is the total decay width of the 
SM Higgs boson with same mass as /i2- We use HDECAY [35] to compute r(ifsM). In 
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the numerical analysis, we find that S' H'H' decay typically has a rate of only 
a few percent, BR{S' H'H') < 5%, and is therefore negligible over most of the 
parameter space. For the mostly-doublet state, on the other hand, Higgs-to-Higgs 
decays can be important in the region niH' < 200 GeV and ms' ~ 50 — 100 GeV. 
Our data set contains six points with rriH' ~ 175 — 200 GeV and ms' ~ 60 — 95 GeV 
which have BR{H' S'S') ~ 3 — 8%. With cos^ > 0.9 these points have small 
doublet-singlet mixing. The singlet here decays 80-85% of the time to bb and about 
9% to r+r^, leading to 46, 262r and 4r final states. We also find one point with 
niH' = 150 GeV, ms' = 64 GeV, cos^ = 0.99 and BR{H' S'S') = 56%. This point 
remains after EWPO contraints. It is also worthwhile to remember that, even when 
Higgs-to-Higgs decays are absent or negligible, the total decay width is modified by 
a factor xi = cos^ 9 (sin^ 6) in the case of H' (S"), relative to the SM Higgs boson. 

Let us finally discuss the discovery potential at the LHC. To this end we use 
the CMS expectations on SM Higgs boson searches presented in [36]. Figure 10.38 
of [36] shows the luminosity needed for a 5a discovery in various standard search 
channels. For H' and 5", this luminosity scales with 1/^^ due to the reduced produc- 
tion cross section, and, where applicable, a factor stemming from the modification 
of the branching ratios into SM particles; for the H': 

BR{H' Xsm) ^r(ij'sM) / , 

BR{HsM ^ Xsm) ~ ^T{Hsm) + TiH' ^ S'S') ' ^ ' ^ 

with ^ = cos^6', and analogously for S' with H' ^ S' and ^ = sin^^^. The resulting 
luminosity needed for a 5cr discovery is shown in figure |1^; the left (right) plot shows 
the lighter (heavier) mass eigenstate. Blue dots represent a mostly-doublet Higgs, 
green dots a mostly-singlet one; points in darker colour are those which survive 
the EWPO constraints. The number and nature of Higgs bosons which are within 
discovery reach with 30 fb^^ of data is shown in figure |18|. 

In summary, we conclude that there is a strong discovery potential for this non- 
standard Higgs sector, if it is the origin of a strong EWPT. This holds in particular for 
the mostly-doublet Higgs, but in many cases also for the mostly-singlet one. There 
are even a couple of points where both H' and S' could be discovered at the LHC with 
30 fb~^ of integrated luminosity. With higher luminosity, the prospects for discover- 
ing both mass eigenstates are quite promising. The alert reader will note, however, 

^Thc branching ratio of S' H'H' is enhanced for cos6' — > 1, but in this case the S' production 
rate goes to zero. Our data set contains one point with ms' = 452 GeV, iriH' ~ 173 GeV, cos 9 — 1 
and BR(5' H'H') = 100%. AU other points have BR{S' H'H') < 5%. 
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Figure 17: Luminosity needed for a 5a discovery of the lighter (left plot) and the heavier 
(right plot) mass eigenstate, extrapolated from CMS results [36]. Blue dots represent 
a mostly-doublet Higgs, green dots a mostly-singlet one; points in darker colour survive 
EWPO constraints. The horizontal dashed lines indicate 30 fb^^ or three years of running 
at low luminosity. 



that figure ^ is limited to 115 GeV < mH',s' ^ 600 GeV. In figure |T8|, we also show 
ms' < 115 GeV, but the S' can be observed in this case. The reason is that the 
experimental analyses in [36] only cover the mass range 115 GeV < rriH < 600 GeV. 
In the majoron model, as well as in other singlet-extensions of the SM, it would 
however be interesting to search also for lighter and heavier states. We hope that 
the present work provides some motivation toward this end. 



5. Conclusions 



We have given an in-depth analysis of the EWPT in the singlet Majoron model at 
the one-loop level, taking account of the LEP constraints on the Higgs boson mass 
and mixing angle. Our broad scan of the model's parameter space indicates that a 
certain level of tuning is needed to get a strong enough transition for electroweak 
baryogenesis; many examples approximately satisfy the relation (|3.71 ) which reduces 
the dimensionality of the parameter space. The zero-temperature one-loop correction 
was shown to play a crucial role in this respect, effectively reducing the negative mass 
squared of the Higgs and hence the critical temperature of the phase transition, to 



25 



600 



500 



400 




300 



200 



100 



Reach with 30 fb"^ 

-I — I — I — I — ! — I — I — r~i — I — I — I — r- 



Reach with 30 fb"' incl. EWPO 



• H' 

• S' 

• both 



1 



J l_l_J ' ' ' ' ' I ' ' I l_L 



600 



500 



400 



300 



200 



100 



J — — 1 — \ — 1 — 1 — \ — 1 — 1 — 1 — r- 
- • • 


T 1 — 1 ] — 1 — 1 — 1 — 1 — ] 1 — 1 ] — 1 — C 

• H' 


- v 


s' ; 




• both : 


- • 




i 




- t • 

- • • 
• 

• 

'I 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 r 1 1 1 1 1 1 1 1 1 r 



100 



200 



300 400 
rriH' [GeV] 



500 



600 



100 



200 



300 400 
rriH' [GeV] 



500 



600 



Figure 18: Discovery reach with 30 fb^^ in the mgi versus niH' plane. The different 
colours encode which state can be detected: blue stands for H' , green for S' , red for both. 
The left plot is without, the right one with EWPO constraints. 



increase the figure of merit for the strength of the transition, Vc/T^. These examples 
required rather large values of the coupling constants, as large as Aj < 3. We also 
identified another population of accepted points with Aj < 1 which has a different 
origin, and exists only for very small singlet masses ms ^ 20 GeV. 

From the technical point of view, our job of identifying cases with a first order 
phase transition was made more difficult by the fact that both fields H and S usually 
evolve during the transition, so it is essential to keep track of both. Although numer- 
ous studies have been done on similar models with a singlet coupling to the Higgs, 
most of these assume that the singlet does not get a VEV, so this complication does 
not arise. Another difference is that a generic model of a real singlet interacting with 
H has many additional couplings which are odd in S, whereas the singlet Majoron 
model is constrained by the U{1) global lepton symmetry, which it spontaneously 
breaks. Due to the reduced number of coupling constants, we were able to make 
an exhaustive search of the parameter space. The singlet Majoron model also has 
the appeal of being theoretically motivated by the seesaw mechanism for neutrino 
masses. In our case, this must be supplemented by the requirement of small Dirac 
Yukawa couplings for the neutrinos, since we take the right-handed neutrinos to be 
lighter than the TeV scale. 

Of course the strong phase transition is only interesting for baryogenesis if there is 
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also a mechanism for producing the baryon asymmetry. Complex Majorana Yukawa 
couplings, which we assumed here to be real for simplicity, could provide the needed 
CP violation. Perhaps CP-violating reflections of the heavy neutrinos at the bubble 
walls could create a lepton asymmetry which would be converted to the baryon 
asymmetry via sphaleron interactions. This is a subject to which we hope to return. 

It would be interesting to extend our study to the generation of gravitational 
waves. Although electroweak baryogenesis and gravity wave generation both need a 
"strong" phase transition, the criteria are different. In particular, relativistic bubble 
walls are favored for producing significant gravity waves. Ref. [33] has recently shown 
that this can be achieved in the related model of a real singlet with more general 
couplings (not respecting any Z2 symmetry) than in the Majoron model. 

Concerning collider phenomenology, we have shown that the LHC has a strong 
discovery potential for this nonstandard Higgs sector, if it is the origin of a strong 
EWPT. This includes the possibility of a SM-like Higgs boson with mass up to 
about 200 GeV, which has a sizeable branching fraction into a pair of light singlets. 
Moreover, with high enough luminosity there are good prospects to discover both 
the H' and the S' states. A dedicated experimental study would be worthwhile to 
cover masses below 115 GeV and above 600 GeV. 
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A. Analytic approximation for thermal potential 

The method of smoothly matching the low- and high-T expansions for the one-loop 
thermal potential was given in ref. [5] . For convenience we repeat the formulas here. 
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The nth order high-T (small M/T) expansion is given by [20] 



Vs,h[n) + - W ~ 64^^ (^^Og \-W)-(^bj 

n J 

:M2ll\^f^ML\ (2/-3)!!C(2^-i) bosons- 

^ 2 l,47r2T2 J (2/)!!(/+l) ' UObOnS, 

1=2 

V„,{n) = -I|gl + Mfl + ^ (log (fl) - 

-^E(3^)'«^(2^'-'-l). felons; 

1=2 

Cb = 3/2 + 2 log 47r- 27b = 5.40762; c/ = - 2 log4 = 2.63503 (A.l) 

respectively for bosons and fcrmions. The corresponding low-T (large M/T) expan- 
sion is [9] 

M(.)^-e-(^)%|:J,Q||±|(r/M, (A.) 

By trial and error, one can find that the low- and high-T expansions can be 
smoothly matched onto each other using the approximation 

H = Q{x, - {M/Tf) y,,b(3) + e((M/r)2 - x,) (M(3) - 5,T') 

Vf = Q{xf - {M/Tf) (4) + e((M/r)2 - xf) m) - 5f T') (A.3) 

where © is the step function with Xi, — 9.47134 and Xf — 5.47281 for bosons and 
fermions, respectively. The small constant shifts of Vi{3) are made so that the func- 
tion as well as its derivatives match at the transition point: Sb — 3.19310 x 10~^ and 
5/ = 4.60156xl0~^. This gives an approximation with a relative error which is less 
than 0.5% for M/T oo, and neghgible for small M/T. 

B. Field dependent masses 

Although it is convenient to express the potential and the field-dependent masses in 
terms of the complex VEV's, it is simpler to compute the masses in the real basis, 
where H° ^ {h + i(/)i)/V2, H+ = (02 + ih)/\/2, 5 = (s + ij)/\/2. We continue to 
express the field-dependent masses in terms of the complex fields, but assuming only 
h and s actually get VEV's: H = h/\/2, S = s/\/2. 
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The diagonal components of the zero-temperature scalar mass matrix, in the real 
basis h, 4>i, s,j, are 







-vI)+\hs\S\' 






-vl)+Xhs\S\' 




= A, {6\S\^ - 


-vl)+Xhs\H\'' 




= A, {2\S2\^ 


-vl)+X,,\H\' 



Within the full 6x6 mass matrix, there is only one off-diagonal entry, 

ml, = ml, = 2X,s\H\\S\ (B.2) 

(Note that we can take H and S to be real here.) Thus one can analytically find all 
the field-dependent mass squared eigenvalues, by diagonalizing the h-s sector. 

For the ring improvement, we must add thermal corrections to the mass squared 
matrix. 




They can be computed by inserting the zero-temperature masses into the high-T 
expansion of the one-loop thermal potential, and reading off the corrections to the 
mass terms (the coefficients of the terms quadratic in H and -S"). These thermal 
masses are the same for each real component {h,4>i or s,j) within the H or S fields, 
respectively. Since Srri^ ^ the thermal mass matrix has to be diagonalized 

independently of the zero-temperature mass matrix. 

The field-dependent masses of the relevant fermions are given by 

m',^y^\Hf, ml^y^lSf (B.5) 

They do not need to be thermally corrected for the ring improvement. In the imag- 
inary time formalism of finite-temperature field theory, where the effective squared 
masses of the Matsubara modes are M^(0, T) + {2'KnTy for bosons and M^(0, T) -\- 
{27: {n + \TY for fermions. Only for the n — Q modes of the bosons can there be 
an infrared divergence due to vanishing M^(0) which would make it important to 
include the perturbative g^T'^ contribution to M^. 
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The only other fields we must consider are the gauge bosons. In the basis of 
Wi, W2, W3, B, the mass matrix is 



fa' 



9'^ 99' 

I „/2 



\ 99' 9" J 



fa' 



\ 



\ 



9" J 



2, longitudinal; 
0, transverse. 



(B.6) 

Only the longitudinal components get a thermal correction at leading order in the 
gauge couplings. 



C. Renormalization constants 



To express the solutions to the renormalization conditions ( |2.5| ) it is convenient to 
define multiplicities gi for the respective fields as: 1 for each real scalar, —12 for the 
top quark, —2 for each of the three right-handed neutrinos, 2 for each transverse 
gauge boson, 1 for longitudinal gauge bosons. It is straightforward to show that 



In/i^ 
A 



i9imi^i}'^mi - I) 
i9^mi^ 



In 



1717 



- 1 



(C.l) 
(C.2) 



evaluated at the minimum of the tree level potential. These conditions ensure that 
the position of this minimum remains unchanged at one loop. 



D. Beta functions 



Defining jSx = IQir'^dX/dlnfi', the beta functions for the largest couplings in the 
singlet Majoron model are [34] 
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(D.l) 




(D.3) 



(D.4) 



(D.5) 



We integrate these starting from the scale /x =100 GeV up to the first Landau 
pole (where any of the running couplings diverge), taking 3 generations of right- 
handed neutrinos. The running of the gauge couplings is neglected in this estimate. 

E. Electroweak precision observables 

To evaluate the impact of constraints on the oblique parameters S, T, U, we follow the 
procedure of references [14,30], defining Ax^ as in eq. (5.4) of [14], and taking points 
with Ax^ > 7.8 to be excluded at 95% c.l. Explicit expressions for T were given in 
those references, but not for S or U. These can be derived from the definitions found 
in eqs. (10.61) of the PDG Review of Particle Properties [37], and the expressions 
for the W and Z self-energies in appendix A of [30] . We find that the contribution 
to S from the Higgs sector is 



-27r5, 



new 



= cos' e [G{H', Z) - ml>Fi{Z, H', Z)] + F2{H' , Z, Z) 



+ 2 [F^iW, H', 0) - Fo{W, H', W) + Fo{Z, H', Z) - Fo{Z, H', 0)] 




(E.l) 



where 



G{a, 6) = - {ml + 



I) 





(E.3) 
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and we define mg = 0. The deviation from tfie standard model prediction is obtained 
by taking 

AS = Sncw - 5'ncw(0 = 0, m^, = m% = ml) (E.4) 
Similarly, U can be inferred from the combination 

[F,{W, H\ 0) - F,{W, H', Z, W)] + F{H\ W, 0) 

+ 2 [Fq{W, H', W) - Fo{W, H', 0)] + F,{H', W, 0) - Fi{H', W, W) 

+ sin^ 9 (^H' ^ S'^ (E.5) 

where F = Fi — F2. AU is computed analogously to (p^. 
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